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GLOSSARY 


Crystallizer The vessel or process unit in which crystal- 
lization occurs. 

Growth The increase in crystal size due to deposition of 
solute on crystal surfaces. 

Magma The mixture of crystals and mother liquor in the 

crystallizer. 

Mode of crystallization The means by which a thermo- 

dynamic driving force for crystallization is created. 

Mother liquor The liquid solution from which crystals 

are formed. 

MSMPR erystallizer A vessel operating in a continuous 
manner in which crystallization occurs and whose con- 
tents are perfectly mixed. As a result of perfect mixing, 
all variables descriptive of the mother liquor and crys- 
tals are constant throughout the vessel and are identical 
to corresponding variables in the product stream leav- 
ing the vessel. 

Nucleation The formation of new crystals. 

Primary nucleation The formation of crystals by mech- 
anisms that do not involve existing crystals of the crys- 


tallizing species; includes both homogeneous and het- 
erogeneous nucleation mechanisms. 

Secondary nucleation The formation of crystals through 
mechanisms involving existing crystals of the crystal- 
lizing species. 

Solubility The equilibrium solute concentration. The di- 
mensions in which solubility is expressed include, but 
are not limited to, mass or mole fraction, mass or mole 
ratio of solute to solvent, and mass or moles of solute 
per unit volume of solvent or solution. 

Supersaturation The difference between existing and 
equilibrium conditions; the quantity represents the 
driving force for crystal nucleation and growth. 


CRYSTALLIZATION PROCESSES addressed in this 
discussion are used in the chemical, petrochemical, phar- 
maceutical, food, metals, agricultural, electronics, and 
other industries. Moreover, the principles of crystalliza- 
tion are important in all circumstances in which a solid 
crystalline phase is produced from a fluid, even when the 
solid is not a product of the process. Much has been done 
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in recent years to improve the understanding of crystal- 
lization, and a large portion of the research on the topic 
has dealt with mechanisms of nucleation and growth. Es- 
pecially important has been elucidation of the effects of 
process variables on the rates at which these phenomena 
occur. Additionally, extensive progress has been achieved 
in modeling both steady-state and dynamic behavior of 
crystallization systems of industrial importance. The pri- 
mary elements of the discussion that follows are the princi- 
ples that influence yield, morphology, and size distribution 
of crystalline products. 


I. OBJECTIVES OF CRYSTALLIZATION 
PROCESSES 


Several examples of objectives that may be satisfied in 
crystallization processes are given in the following dis- 
cussion. Soda ash (sodium carbonate) is recovered from 
brine by contacting the brine with carbon dioxide that re- 
acts with sodium carbonate to form sodium bicarbonate. 
Sodium bicarbonate, which has a lower solubility than 
sodium carbonate, crystallizes as it is formed. The pri- 
mary objective of the crystallizers used in this process 
is separation of a high percentage of sodium bicarbonate 
from the brine in a form that facilitates segregation of the 
crystals from the mother liquor. The economics of crystal 
separation from the mother liquor are affected primarily 
by the variables that control the flow of liquid through the 
cake of crystals formed on a filter or in a centrifuge. For 
example, the flow rate of liquid through a filter cake de- 
pends on the available pressure drop across a filter, liquid 
viscosity, and the size distribution of crystals collected on 
the filter. With a fixed available pressure drop and defined 
liquid properties, the crystal size distribution controls filter 
throughput and, concomitantly, the production rate from 
the process. 

Crystallization can be used to remove solvent from a 
liquid solution. For example, concentration of fruit juice 
requires the separation of solvent (water) from the natural 
juice. The common procedure is evaporation, but the de- 
rived juices may lose flavor components or undergo ther- 
mal degradation during the evaporative process. In freeze 
concentration, the solvent is crystallized (frozen) in rela- 
tively pure form to leave behind a solution with a higher 
solute concentration than the original mixture. Significant 
advantages in product taste have been observed in the ap- 
plication of this process to concentrations of various types 
of fruit juice. 

The elimination of small amounts of an impurity from a 
product species may be an objective of crystallization. In 
such instances, a multistep crystallization—redissolution— 
recrystallization process may be required to produce a 
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product that meets purity specifications. For example, in 
the manufacture of the amino acid L-isoleucine, the prod- 
uct is first recovered in acid form, redissolved, neutralized, 
and then recrystallized in order to exclude the impurity 
L-leucine and other amino acids from the product. 

A simple change in physical properties also can be 
achieved by crystallization. In the process of making soda 
ash, referred to earlier, the sodium bicarbonate crystals are 
subjected to heat that causes the release of carbon dioxide 
and produces low-density sodium carbonate crystals. The 
density of these crystals is incompatible with their use in 
glass manufacture, but a more acceptable crystal can be 
obtained by contacting the sodium carbonate crystals with 
water to form crystalline sodium carbonate monohydrate. 
Drying the resulting crystals removes the water of hydra- 
tion and produces a dense product that is acceptable for 
glass manufacture. 

Separation of a chemical species from a mixture of sim- 
ilar compounds can be achieved by crystallization. The 
mode of crystallization may fall in the realm of what 
is known as melt crystallization. In such processes, the 
mother liquor largely is comprised of the melt of the crys- 
tallizing species, and, subsequent to its crystallization, 
crystals formed from the mother liquor are remelted to 
produce the product from the crystallizer. Para(p)-xylene 
can be crystallized from a mixture that includes ortho 
and meta isomers in a vertical column that causes crys- 
tals and mother liquor to move countercurrently. Heat is 
added at the bottom of the column to melt the p-xylene 
crystals; a portion of the melt is removed from the crys- 
tallizer as product and the remainder flows up the col- 
umn to contact the downward-flowing crystals. Effluent 
mother liquor, consisting almost entirely of the ortho and 
meta isomers of xylene, is removed from the top of the 
column. 

Production of a consumer product in a form suitable 
for use and acceptable to the consumer also may be an ob- 
jective of a crystallization process. For example, sucrose 
(sugar) can be crystallized in various forms. However, dif- 
ferent cultures are accustomed to using sugar that has a 
particular appearance and, unless the commodity has that 
appearance, the consumer may consider the sugar to be 
unacceptable. 

In all these processes, there are commmon needs: to 
form crystals, to cause them to grow, and to separate the 
crystals from the residual liquid. While conceptually sim- 
ple, the operation of a process that utilizes crystallization 
can be very complex. The reasons for such complexity in- 
volve the interaction of the common needs and process re- 
quirements on product yield, purity, and, uniquely, crystal 
morphology and size distribution. In the following discus- 
sion, the interactions will be discussed and general princi- 
ples affecting crystallizer operation will be outlined. More 
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extensive discussion of the subject matter can be found in 
the bibliography at the end of the chapter. 


ll. EQUILIBRIUM AND MASS AND 
ENERGY BALANCES 


A. Solid—Liquid Equilibrium 


The solubility of a chemical species in a solvent refers 
to the amount of solute that can be dissolved at constant 
temperature, pressure, and solvent composition (including 
the presence of other solutes). In other words, it is the 
concentration of the solute in the solvent at equilibrium. 

As with all multiphase systems, the Gibbs phase rule 
provides a useful tool for determining the number of inten- 
sive variables (ones that do not depend on system mass) 
that can be fixed independently: 


Np = Ne — Np + 2 (1) 


Nor is the number of degrees of freedom, Ne is the num- 
ber of components, and N, is the number of phases in the 
system. The number of degrees of freedom represents the 
number of independent variables that must be specified 
in order to fix the condition of the system. For example, 
the Gibbs phase rule specifies that a two-component, two- 
phase system has two degrees of freedom. If temperature 
and pressure are selected as the specified variables, then 
all other intensive variables—in particular, the composi- 
tion of each of the two phases—are fixed, and solubility 
diagrams of the type shown for a hypothetical mixture of 
R and S in Fig. 1 can be constructed. 

Several features of the hypothetical system described 
in Fig. | illustrate the selection of crystallizer operating 
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FIGURE 1 Hypothetical solubility diagram of eutectic-forming 
system. 
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conditions and the limitations placed on the operation by 
the system properties. The curves AB and BC represent 
solution compositions that are in equilibrium with solids 
whose compositions are given by the lines AD and CE, 
respectively. If AD and CE are vertical and coincident 
with the left and right extremes, the crystals are pure S and 
R, respectively. Crystallization from any solution whose 
equilibrium composition is to the left of a vertical line 
through point B will produce crystals of pure S, while 
solutions with an equilibrium composition to the right of 
the line will produce crystals of pure R. A solution whose 
composition falls on the line through B will produce a 
mixture of crystals of R and S. 

Now suppose a saturated solution at temperature T; is 
fed to a crystallizer operating at temperature Tz. Since it is 
saturated, the feed has a mole fraction of R equal to xp. The 
maximum production rate of crystals occurs when the so- 
lution leaving the crystallizer is saturated, meaning that the 
crystal production rate, mproa, depends on the value of 7): 


Mprod = MFXF — MLXL (2) 


where mp is the feed rate to the crystallizer and m, is the 
solution flow rate leaving the crystallizer. Note that the 
lower limit on T is given by the eutectic point, and that 
attempts to operate the crystallizer at a temperature other 
than the eutectic value will result in a mixture of crystals 
of R and S. 

When certain solutes crystallize from aqueous solu- 
tions, the crystals are hydrated salts, which means that 
the crystals contain water and solute in a specific stoichio- 
metric ratio. The water in such instances is referred to as 
water of hydration, and the number of water molecules 
associated with each solute molecule may vary with the 
crystallization temperature. 

Potassium sulfate provides an example of such behav- 
ior. When it crystallizes from an aqueous solution above 
40°C, the crystals are anhydrous K2SO,, while below 
40°C each molecule of K2SO, that crystallizes has 10 
molecules of water associated with it. The hydrated salt, 
K2SO,-10H20(s), is called potassium sulfate decahydrate. 
Another solute that forms hydrated salts is magnesium 
sulfate, which can incorporate differing amounts of water 
depending upon the temperature at which crystallization 
occurs (see Table I). 

The solubility diagrams of several species are shown 
in Fig. 2, and these illustrate the importance of solubility 
behavior in the selection of the mode of crystallization. 
For example, consider the differences between potassium 
nitrate and sodium chloride: The solubility of potassium 
nitrate is strongly influenced by the system temperature, 
whereas the opposite is true for sodium chloride. As a 
consequence, (1) a high yield of potassium nitrate crys- 
tals can be obtained by cooling a saturated feed solution, 
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TABLE! Water of Hydration for MgSO, 

Form Name wt% MgSO4 Conditions 
MgSO4 Anhydrous magnesium sulfate 0.0 >100°C 
MgSO4-H20 Magnesium sulfate monohydrate 87.0 67 to 100°C 
MgS04-6 H20 Magnesium sulfate hexahydrate 52.7 48 to 67°C 
MgS04-7 H2O Magnesium sulfate heptahydrate 48.8 2 to 48°C 
MgS04-12 H20 Magnesium sulfate dodecahydrate 35.8 —4 to 2°C 


but (2) cooling a saturated sodium chloride solution ac- 
complishes little crystallization, and vaporization of water 
is required to increase the yield. 

The effect of water of hydration on solubility can be seen 
in Fig. 2. Note, for example, that sodium sulfate has two 
forms in the temperature range of the solubility diagram: 
sodium sulfate decahydrate (Na2S04-10H20), which is 
known as Glauber’s salt, and anhydrous sodium sulfate. 
Since a transition from Glauber’s salt to anhydrous sodium 
sulfate occurs at approximately 34°C, crystals recovered 
from a crystallizer operating above about 34°C will be 
anhydrous, but those from a crystallizer operating below 
this temperature will contain 10 waters of hydration. Also 
observe the effect of water of hydration on solubility char- 
acteristics; clearly, cooling crystallization could be used to 
recover significant yields of Glauber’s salt but evaporative 
crystallization would be required to obtain high yields of 
the anhydrous salt. 

Mixtures of multiple solutes in a single solvent are en- 
countered in a number of processes—for example, in the 
recovery of various chemicals from ores or brines. Expres- 
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sion of the complex solubility behavior in such systems 
by graphical means usually is limited to systems of two 
solutes. The interaction of added solutes on solubility is 
illustrated by the plot of equilibrium behavior for potas- 
sium nitrate—sodium nitrate—water in Fig. 3. As before, 
the curves in the diagram trace solution compositions that 
are in equilibrium with solid solutes. Points A, D, G, and 
J are based on the solubilities of pure potassium nitrate, 
while points C, F, I, and L are based on solubilities of 
pure sodium nitrate. Curves AB, DE, GH, and JK repre- 
sent compositions of solutions in equilibrium with solid 
potassium nitrate at 30, 50, 70, and 100°C, respectively. 
Curves BC, EF, HI, and KL represent compositions of so- 
lutions in equilibrium with solid sodium nitrate. Should 
the solution condition, including temperature, correspond 
to points B, E, H, K or any condition on the curve con- 
necting these points, crystals of both solutes would be 
formed by cooling. 

A second type of solubility behavior is exhibited by 
mixtures that form solid solutions. Consider, for exam- 
ple, a hypothetical system containing R and S whose 
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FIGURE 2 Solubility diagram for several common substances. 
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KNO, NaNO, 


FIGURE 3 Solubility diagram of KNO3 and NaNO; mixtures in 
water. 


equilibrium behavior is described in Fig. 4. The phase en- 
velope is drawn based on the compositions of coexisting 
liquid and solid phases at equilibrium. The pure compo- 
nent R has a melting point at pressure P equal to T) while 
the melting point of pure S is Tı. The system behavior can 
best be described by the following example: Consider a 
mixture of R and S at temperature Tą and having a mass 
fraction of R equal to zm. From the phase diagram, the 
mixture is a liquid. As the liquid is cooled, a solid phase 
forms when the temperature reaches Tg and the system is 
allowed to come to equilibrium; the solid-phase composi- 
tion corresponds to a mass fraction of R equal to xg. On 
cooling the liquid further, the ratio of solid to liquid in- 


Mass fraction R 


FIGURE 4 Hypothetical solubility diagram of mixture without a 
eutectic at constant pressure: x, solid; y, liquid; z, combined. 
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creases and at Tc the mass fraction of R in the liquid is yc 
and in the solid it is xc. At Tp the liquid phase disappears, 
leaving a solid with a mass fraction of R equal to zy. 

Systems that exhibit behavior of the type illustrated 
in Fig. 4 cannot be purified in a single crystallization 
stage. They represent situations in which multiple stages 
or continuous-contacting devices may be useful. The prin- 
ciples of such operations are analogous to those of other 
countercurrent contacting operations—for example, dis- 
tillation, absorption, and extraction. 

Variables other than temperature and the presence of 
other solutes can influence solubility. For example, the ef- 
fect of a nonsolvent on solubility sometimes is used to 
bring about recovery of a solute. Figure 5 shows the sol- 
ubility of L-serine in aqueous solutions containing vary- 
ing amounts of methanol. Note that increasing methanol 
content reduces the solubility by more than an order of 
magnitude, and this characteristic can be used to obtain a 
high yield in the recovery of L-serine. 

There is increasing interest in the crystallization of so- 
lutes from supercritical-fluid solvents. In such instances, 
solubilities often are correlated by an equation of state. 
Such concepts are beyond the scope of the current discus- 
sion but are presented elsewhere in the encyclopedia. 

Although this discussion provides insight to the types 
of solubility behavior that can be exhibited by various sys- 
tems, itis by no means a complete survey of the topic. Ex- 
tensive solubility data and descriptions of more complex 
equilibrium behavior can be found in the literature. Pub- 
lished data usually consist of the influence of temperature 
on the solubility of a pure solute in a pure solvent; seldom 
are effects of other solutes, co-solvents, or pH considered. 
As a consequence, solubility data on a system of interest 
should be measured experimentally, and the solutions used 
in the experiments should be as similar as possible to those 
expected in the process. Even if a crystallizer has been de- 
signed and the process is operational, obtaining solubility 
data using mother liquor drawn from the crystallizer or 
a product stream would be wise. Moreover, the solubility 
should be checked periodically to see if it has changed due 
to changes in the upstream operations or raw materials. 

There have been advances in the techniques by which 
solid—liquid equilibria can be correlated and, in some 
cases, predicted. These are described in references on 
phase-equilibrium thermodynamics. 


B. Mass and Energy Balances 


Illustrating the formulation of mass and energy bal- 
ances is simplified by restricting the analysis to sys- 
tems whose crystal growth kinetics are sufficiently fast 
to utilize essentially all of the supersaturation provided 
by the crystallizer; in other words, the product solution 
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FIGURE 5 Effect of methanol on solubility of L-serine. 


is assumed to be saturated. Under such conditions (re- 
ferred to in the crystallization literature as Class II or fast- 
growth behavior), the solute concentration in the mother 
liquor can be assigned a value corresponding to satura- 
tion. Should the supersaturation in the mother liquor be 
so great as to affect the solute balance, the operation is 
said to follow Class I or slow-growth behavior. In Class I 
behavior, the operating conditions affect the rate at which 
solute is crystallized, and an expression coupling the rate 
of growth to a solute balance must be used to describe 
the system. Such treatment will be considered beyond the 
scope of this discussion. 

The solution of mass and energy balances requires sol- 
ubility and enthalpy data on the system of interest. Various 
methods of presenting solubility data were given earlier, 
and the use of solubilities to estimate crystal production 
rates from a cooling crystallizer was demonstrated by the 
discussion of Eq. (2). Subsequent to determining the yield, 
the rate at which heat must be removed from the crystal- 
lizer can be calculated from an energy balance: 


mcÊc+m Ê, —mpHp = Q (3) 


where mp, mc, and m, are feed rate, crystal production 
rate, and mother liquor flow rate, respectively; H is spe- 
cific enthalpy of the stream corresponding to the subscript; 
and Q is the required rate of heat transfer to the crystallizer. 
AS mp, Mc, and m, are known or can be calculated from a 
simple mass balance, determination of Q requires estima- 
tion of specific enthalpies. These are most conveniently 
obtained from enthalpy-composition diagrams, which are 
available in the general literature for a number of sub- 
stances. 


If specific enthalpies are unavailable, they can be es- 
timated based on defined reference states for both solute 
and solvent. Often the most convenient reference states are 
crystalline solute and pure solvent at an arbitrarily chosen 
reference temperature. The reference temperature selected 
usually corresponds to that at which the heat of crystal- 
lization, AH c, Of the solute is known. (The heat of crystal- 
lization is approximately equal to the negative of the heat 
of solution.) For example, if the heat of crystallization is 
known at Tef, then reasonable reference conditions would 
be the solute as a solid and the solvent as a liquid, both at 
Tet. The specific enthalpies could be estimated then as: 


Ag = xpAH, + Cp.(T =, Tief) (4) 
He = Cpc(T — Tret) (5) 
Ay = xi AH. + Cp (T = Tref) (6) 


where xp and x; are the mass fractions of solute in the 
feed and mother liquor, respectively. All that is required 
now to determine the required rate of heat transfer is the 
indicated heat capacities, which can be estimated based 
on system composition or measured experimentally. 

Now suppose some of the solvent is evaporated in the 
crystallizer. Independent balances can be written on total 
and solute masses: 


mp = My + m, + mc (7) 
XFMF = XML + XcCMc (8) 


Assuming that the streams leaving the crystallizer are in 
equilibrium, there is a relationship between the temper- 
ature (or pressure) at which the operation is conducted 
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and x, and xc. In addition, an energy balance must be 
satisfied: 


mpHp + Q =myHy +m AL +mcHc (9) 


The specific enthalpies in the above equation can be deter- 
mined as described earlier, provided the temperatures of 
the product streams are known. Evaporative cooling crys- 
tallizers (described more completely in Section V) operate 
at reduced pressure and may be considered adiabatic. In 
such circumstances, Eq. (9) is modified by setting Q = 0. 
As with many problems involving equilibrium relation- 
ships and mass and energy balances, trial-and-error com- 
putations are often involved in solving Eqs. (7) through (9). 


lil. NUCLEATION AND GROWTH KINETICS 


The kinetics of crystallization have constituent phenom- 
ena in crystal nucleation and growth. The rates at which 
these occur are dependent on driving forces (usually ex- 
pressed as supersaturation), physical properties, and pro- 
cess variables, but relationships between these quantities 
and crystallization kinetics often are difficult to express 
quantitatively. As a result, empirical or qualitative links 
between a process variable and crystallization kinetics are 
useful in providing guidance in crystallizer design and op- 
eration and in developing strategies for altering the prop- 
erties of crystalline products. 

Nucleation and growth can occur simultaneously in a 
supersaturated environment, and the relative rates at which 
these occur are primary determinants of the characteris- 
tics of the crystal size distribution; one way of influencing 
product size distributions is through the control of vari- 
ables such as supersaturation, temperature, and mixing 
characteristics. Obviously, those factors that increase nu- 
cleation rates relative to growth rates lead to a crystal size 
distribution consisting of smaller crystals. In the discus- 
sion that follows, an emphasis will be given to the general 
effects of process variables on nucleation and growth, but 
the present understanding of these phenomena does not 
allow quantitative a priori prediction of the rates at which 
they occur. 


A. Supersaturation 


Supersaturation is the thermodynamic driving force for 
both crystal nucleation and growth; and therefore, it is 
the key variable in setting the mechanisms and rates by 
which these processes occur. It is defined rigorously as the 
deviation of the system from thermodynamic equilibrium 
and is quantified in terms of chemical potential, 


Ami = mi — u? = RT In & (10) 
a 


i 
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where w; is the chemical potential of solute i at the existing 
conditions of the system, u¥ is the chemical potential of the 
solute equilibrated at the system conditions, and a; and a} 
are activities of the solute at the system conditions and at 
equilibrium, respectively. Less abstract definitions involv- 
ing measurable system quantities are often used to approx- 
imate supersaturation; these involve either temperature or 
concentration (mass or moles of solute per unit volume or 
mass of solution or solvent) or mass or mole fraction of 
solute. Recommendations have been made that it is best to 
express concentration in terms of moles of solute per unit 
mass of solvent. For systems that form hydrates, the so- 
lute should include the water of hydration, and that water 
should be deducted from the mass of solvent. 

Consider, for example, a system at temperature T with 
a solute concentration C, and define the equilibrium tem- 
perature of a solution having a concentration C as T* and 
the equilibrium concentration of a solution at T as C*. 
These quantities may be used to define the following ap- 
proximate expressions of supersaturation: 


1. The difference between the solute concentration and 
the concentration at equilibrium, AC; = C; — C; 

2. For a solute whose solubility in a solvent increases 
with temperature, the difference between the 
temperature at equilibrium and the system 
temperature, AT = T* — T 

3. the supersaturation ratio, which is the ratio of the 
solute concentration and the equilibrium 
concentration, S; = C; / C¥ 

4. The ratio of the difference between the solute 
concentration and the equilibrium concentration to 
the equilibrium concentration, o; = (C; — C¥)/ C¥ = 
S; — 1, which is known as relative supersaturation. 


Any of the above definitions of supersaturation can be used 
over a moderate range of system conditions, but as outlined 
in the following paragraph, the only rigorous expression 
is given by Eq. (10). 

The definitions of supersaturation ratio and relative su- 
persaturation can be extended to any of the other variables 
used in the definition of supersaturation. For example, 
defining Sa; =a; /a; gives: 


Ali viCi 
RE =InS,, =In eG (11) 
Therefore, for ideal solutions or for y; © y;*, 
ALi Ci 
Hi wm = hns; (12) 
RT C 
Furthermore, for low supersaturations (say, $; < 1.1), 
Ali 
aE gee ee (13) 


RT 
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The simplicity of Eq. (13) results in the use of relative 
supersaturation in most empirical expressions for nucle- 
ation and growth kinetics. While beguilingly simple, and 
correct in limiting cases, great care should be taken in ex- 
tending such expressions beyond conditions for which the 
correlations were developed. 

For ionic solutes, a;=a}ł, which leads to Sy, = 
(a+ /a%)” and 


ALi Vit Cj 
aT = vin Sy, = vin CF (14) 
Again, for yi © y» 
A i Ci 
H vln — = vln $ (15) 
RT CF 


B. Primary Nucleation 


The term primary nucleation is used to describe both ho- 
mogeneous and heterogeneous nucleation mechanisms in 
which solute crystals play no role in the formation of new 
crystals. Primary nucleation mechanisms involve the for- 
mation of crystals through a process in which constituent 
crystal units are stochastically combined. Both homo- 
geneous and heterogeneous nucleation require relatively 
high supersaturations, and they exhibit a high-order de- 
pendence on supersaturation. As will be shown shortly, 
the high-order dependence has a profound influence on 
the character of crystallization processes in which primary 
nucleation is the dominant means of crystal formation. 

The classical theoretical treatment of primary nucle- 
ation that produces a spherical nucleus results in the ex- 
pression: 


o 167 €),,¢0° 
B° = Aexp ; 
3k3T3[In(o + 1)]? 
REEN 16r es pu’ (16) 
Eo 3K3 T 3.02 


where k is the Boltzmann constant, €surf is the interfacial 
surface energy per unit area, v is molar volume of the 
crystallized solute, and A is a constant. 

The theory shows that the most important variables af- 
fecting the rates at which primary nucleation occur are 
interfacial energy €surr, temperature T, and supersatura- 
tion o. The high-order dependence of nucleation rate on 
these three variables, especially supersaturation, is impor- 
tant because, as shown by an examination of Eq. (16), a 
small change in any of the three variables could produce an 
enormous change in nucleation rate. Such behavior gives 
rise to the often observed phenomenon of having a clear 
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liquor transformed to a slurry of very fine crystals with 
only a slight increase in supersaturation, for example by 
decreasing the solution temperature. 

The effect of exogenous solid matter (as in heteroge- 
neous nucleation) in the supersaturated solution is equiv- 
alent to that of a catalyst in a reactive mixture. Namely, 
it is to reduce the energy barrier to the formation of a 
new phase. In effect, the solid matter reduces the interfa- 
cial energy €sy,¢ by what may amount to several orders of 
magnitude. 

The classical nucleation theory embodied in Eq. (16) 
has a number of assumptions and physical properties 
that cannot be estimated accurately. Accordingly, empir- 
ical power-law relationships involving the concept of a 
metastable limit have been used to model primary nucle- 
ation kinetics: 


B° = kyo” (17) 


max 
where ky and n are parameters fit to data and Omax is 
the supersaturation at which nuclei are observed when 
the system is subjected to a specific protocol. Although 
Eq. (17) is based on empiricism, it is consistent with the 
more fundamental Eq. (16). 


C. Secondary Nucleation 


Secondary nucleation is the formation of new crystals 
through mechanisms involving existing solute crystals; in 
other words, crystals of the solute must be present for 
secondary nucleation to occur. Several features of sec- 
ondary nucleation make it important in the operation of 
industrial crystallizers: First, continuous crystallizers and 
seeded batch crystallizers have crystals in the magma that 
can participate in secondary nucleation mechanisms. Sec- 
ond, the requirements for the mechanisms of secondary 
nucleation to be operative are fulfilled easily in most in- 
dustrial crystallizers. Finally, many crystallizers are oper- 
ated in a low supersaturation regime so as to maximize 
yield, and at such supersaturations the growth of crystals 
is more likely to produce desired morphologies and high 
purity; these low supersaturations can support secondary 
nucleation but not primary nucleation. 


1. Mechanisms 


Secondary nucleation can occur through several mecha- 
nisms, including initial breeding, contact nucleation (also 
known as collision breeding), and shear breeding. Al- 
though a universal expression for the kinetics of secondary 
nucleation does not exist, a working relationship often can 
be obtained by correlating operating data from a crystal- 
lizer with a semi-empirical expression. Guidance as to the 
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form of the expression and the variables that it should in- 
clude can be obtained by understanding the various mech- 
anisms of secondary nucleation. 

Initial breeding results from immersion of seed crys- 
tals in a supersaturated solution, and it is thought to be 
caused by dislodging extremely small crystals that were 
formed on the surface of larger crystals during drying. Al- 
though this mechanism is unimportant in continuous and 
unseeded batch crystallization, it can have a significant 
impact on the operation of seeded batch crystallizers. The 
number of crystals formed by initial breeding, has been 
found to be proportional to the surface area of crystals 
used to seed a batch crystallizer. Characteristics of the 
resulting distribution are affected strongly by the growth 
kinetics of nuclei resulting from initial breeding, and the 
phenomenon of growth-rate dispersion (which will be dis- 
cussed later) can lead to erroneous conclusions regarding 
the nucleation kinetics. 

Shear breeding results when supersaturated solution 
flows by acrystal surface and carries with it crystal precur- 
sors believed formed in the region of the growing crystal 
surface. High supersaturation is required for shear breed- 
ing to produce significant numbers of nuclei. 

Contact nucleation in industrial processes results from 
collisions of crystals with the impeller used for circulation 
of the magma or with other crystallizer internals such as 
baffles, pipe and crystallizer walls, and even other crystals. 
Careful experimental studies have shown that the num- 
ber of crystals produced by collisions between crystals 
and these objects depends upon the collision energy, su- 
persaturation at impact, supersaturation at which crystals 
mature, material of the impacting object, area and angle 
of impact, and system temperature. The collision energy 
for contact nucleation is small and does not necessarily 
result in the macroscopic degradation or attrition of the 
contacted crystal. 

Nucleation from collisions between crystals in the cir- 
culating magma and the rotor in a circulation pump or an 
agitator usually dominate nucleation resulting from other 
collisions. The operating variables in systems of this type 
can be manipulated to some extent, thereby modifying 
nucleation rates and the concomitant crystal size distri- 
bution. For example, internal classification can be used to 
keep larger crystals away from energetic collisions with an 
impeller, but doing so may create other problems with sta- 
bility of the crystal size distribution. The rotational speed 
of an impeller can be changed if there are appropriate con- 
trols on the pump or agitator. Caution must be exercised, 
however, for a reduction in circulation velocity can reduce 
heat-transfer coefficients and increase fouling (encrusta- 
tion) on heat-transfer surfaces. Moreover, the crystals in 
the magma must be kept suspended or crystal morphol- 
ogy and growth rates could be affected adversely. Impact 
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energy may have a high-order dependence on rotational 
speed and, if that is the case, modest changes in this vari- 
able could alter nucleation rates substantially. The fraction 
of the impact energy transmitted from an impeller to the 
crystal can be manipulated by changing the material of 
construction of the impeller. The influence of using soft 
materials to coat impellers or crystallizer internals may 
vary from one crystalline system to another; those systems 
in which the crystal face is soft may be more susceptible 
to nucleation rate changes than those crystalline systems 
where the face is hard. 

Supersaturation has been observed to affect contact nu- 
cleation, but the mechanism by which this occurs is not 
clear. There are data that infer a direct relationship be- 
tween contact nucleation and crystal growth; these data 
showed that the number of nuclei produced by an impact 
was proportional to the linear growth rate of the impacted 
face. This could indicate that the effect of supersaturation 
is to alter growth rates and, concomitantly, the charac- 
teristics of the impacted crystal faces; alternatively, what 
appears to be a mechanistic relationship actually could be 
a result of both nucleation and growth depending upon 
supersaturation. 

Another theory that could account for the effect of su- 
persaturation on contact nucleation is based on the view 
that nuclei formed cover a range of sizes that includes the 
critical nucleus. Since only the nuclei larger than the crit- 
ical nucleus are stable, the relationship of the size of the 
critical nucleus to supersaturation reflects the dependence 
of contact nucleation on supersaturation. This concept, 
which has been referred to as a survival theory, seems to 
have been refuted by measurements of the sizes of crystals 
formed by collisions. These sizes are much larger than the 
critical nucleus, and the survival theory would have little 
influence on the number of nuclei that survive. 

Evidence of the formation of polymolecular clusters in 
supersaturated solutions may provide a mechanistic inter- 
pretation of the effect of supersaturation on contact nu- 
cleation kinetics. These clusters may participate in nucle- 
ation, although the mechanism by which this would occur 
is not clear. One model that has been proposed, however, 
calls for the formation of a semi-ordered region consisting 
of molecular clusters awaiting incorporation into the crys- 
tal lattice. Collisions or fluid shear of the region containing 
high cluster concentrations could then result in these clus- 
ters serving as secondary nuclei. In such a model, the vari- 
ables that influence formation and diffusion of the clusters 
also influence crystal growth rates and nucleation. 


2. Kinetic Expressions 


Irrespective of the actual mechanisms by which con- 
tact nucleation occurs, empirical power-law expressions 
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provide a useful means of correlating nucleation kinetics 
and using the resulting correlations in process analysis and 
control. The correlations generally take the form: 


B° = kyo’ MÅ NÝ (18) 


where ky, i, j, k are positive parameters obtained from 
data correlation, Mr is the magma density (mass of solids 
per unit volume of slurry or solvent in the magma), and N 
is the rotational velocity of the impeller or pump rotor. For 
convenience, either crystal growth rate or mean residence 
time, both of which are directly related to supersaturation, 
may be substituted for o in Eq. (18). 

If primary nucleation dominates the process, i tends 
to larger values (say greater than 3), j and k approach 
zero, and Eq. (18) approaches Eq. (17). Should crystal- 
impeller and/or crystal—crystallizer impacts dominate, j 
approaches 1; on the other hand, if crystal—crystal contacts 
dominate, j approaches 2. 

The ease with which nuclei can be produced by con- 
tact nucleation is a clear indication that this mechanism 
is dominant in many industrial crystallization operations. 
Research on this nucleation mechanism is continuing with 
the objective of building an understanding of the phe- 
nomenon that will allow its successful inclusion in models 
describing commercial systems. 


D. Fundamentals of Crystal Growth 


Crystal growth rates may be expressed as (1) the linear 
advance rate of an individual crystal face, (2) the change 
in a characteristic dimension of a crystal, or (3) the rate of 
change in mass of acrystal or population of crystals. These 
different expressions are related through crystal geometry; 
it is often convenient to use the method of measurement 
as the basis of the growth rate expression or, in certain 
instances, the method used to analyze a crystallization 
process will require that growth rate be defined in a spe- 
cific way. For example, the use of a population balance 
to describe crystal size distribution requires that growth 
rate be defined as the rate of change of a characteristic 
dimension. 

Single-crystal growth kinetics involve the advance rate 
of an individual crystal face normal to itself or the rate of 
change in crystal size associated with exposure to a super- 
saturated solution. The advance rate of a single crystal face 
can be quantified by observation of the face through a cali- 
brated eyepiece of an optical microscope, which allows ex- 
amination of the structure of the advancing crystal face and 
isolation of surface-reaction kinetics from mass-transfer 
kinetics (these phenomena will be discussed later). An 
additional advantage of single-crystal systems is that it is 
possible to examine crystal growth kinetics without inter- 
ference from competing processes such as nucleation. 
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Multicrystal-magma studies usually involve examina- 
tion of the rate of change of a characteristic crystal di- 
mension or the rate of increase in the mass of crystals in 
a magma. The characteristic dimension in such analyses 
depends upon the method used in the determination of 
crystal size; for example, the second largest dimension is 
measured by sieve analyses, while an equivalent spherical 
diameter is determined by both electronic zone sensing 
and laser light scattering instruments. A relationship be- 
tween these two measured dimensions and between the 
measured quantities and the actual crystal dimensions can 
be derived from appropriate shape factors. Volume and 
area shape factors are defined by the equations: 


Verys = kaL? and Aerys = kral” (19) 


where Verys and acrys are volume and area of a crystal, 
kyo and karea are Volume and area shape factors, and L 
is the characteristic dimension of the crystal. Suppose an 
equivalent spherical diameter L phere is obtained from an 
electronic zone-sensing instrument, and the actual dimen- 
sions of the crystal are to be calculated. Assume for the 
sake of this example that the crystals have a cubic shape. 
Let Lecube be the edge length of the crystal and ce and 
kebe be the volume shape factors for a sphere and a cube, 
respectively. Since the volume of the crystal is the same, 
regardless of the arbitrarily defined characteristic dimen- 
sion, 


sphere y 3 = cube 7 3 (20) 


Verys = Kol sphere vol ““cube 


Since kP? is 2/6 and k®®™ is 1.0, the numerical rela- 


tionship between Leube and Lsphere is given by: 


sphere 1/3 1/3 
Kol T 
Loube = pevbe L sphere = A L sphere (21) 


vol 


The rate of change of a crystal mass dMerys/dt can be 
related to the rate of change in the crystal characteristic 
dimension (dL /dt = G) by the equation: 


dMeys d(pkyoiL*) 2 dL 
Z= = 3pky L* | — 22 
di dt AECA A a 
where p is crystal density. Since karea = Acrys/ EZ 
dMexys 
a = 30(Kkvoi/ Karea)AcrysG (23) 


At least two resistances contribute to the kinetics of 
crystal growth. These resistances apply to (1) integration 
of the crystalline unit (e.g., solute molecules) into the crys- 
tal surface (1.e., lattice), and (2) molecular diffusion or bulk 
transport of the unit from the surrounding solution to the 
crystal surface. As aspects of molecular diffusion and mass 
transfer are covered elsewhere, the current discussion will 
focus only on surface incorporation. 
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1. Mechanisms 


Among the many models that have been proposed to de- 
scribe surface-reaction kinetics are those that assume crys- 
tals grow by layers and others that consider growth to oc- 
cur by the movement of a continuous step. Each physical 
model results in a specific relationship between growth 
rate and supersaturation and, although none can predict 
growth kinetics a priori, insights regarding the effects of 
process variables on growth can be obtained. Because of 
the extensive literature on the subject, only the key aspects 
of the physical models and (in one case) the resulting re- 
lationship between growth and supersaturation predicted 
by each theory will be discussed here. 

The model used to describe the growth of crystals by 
layers is based on a two-step, birth-and-spread mecha- 
nism. In one of the steps (birth) a two-dimensional nu- 
cleus is formed on the crystal surface, and in the sec- 
ond step (spread) the two-dimensional nucleus grows to 
cover the crystal surface. When one or the other of the 
steps is controlling growth rates, simplifications of the 
more complicated dependence of growth rate on super- 
saturation can be developed to give what are known as 
the mononuclear two-dimensional nucleation theory and 
the polynuclear two-dimensional nucleation theory. 
In the mononuclear two-dimensional nucleation theory, 
surface nucleation occurs at a finite rate while the spread- 
ing across the surface occurs at an infinite rate. The re- 
verse is true for the polynuclear two-dimensional nucle- 
ation theory. Theoretical relationships have been derived 
between growth rate and supersaturation for each of these 
conditions but are considered beyond the scope of this 
discussion. 

The screw-dislocation theory (sometimes referred to as 
the BCF theory because of its development by Burton, 
Cabrera, and Frank) is based on a mechanism of continu- 
ous movement in a spiral or screw of a step or ledge on the 
crystal surface. The theory shows that the dependence of 
growth rate on supersaturation can vary from a parabolic 
relationship at low supersaturations to a linear relationship 
at high supersaturations. In the BCF theory, growth rate is 


given by: 
co? b 
G = kg{ — } tanh | — (24) 
b Eo 


where e€ is screw dislocation activity and b is a system- 
dependent quantity that is inversely proportional to tem- 
perature. It can be shown that the dependence of growth 
rate on supersaturation is linear if the ratio b/e€o is large, 
but the dependence becomes parabolic as the ratio be- 
comes small. It is possible, then, to observe variations in 
the dependence of growth rate on supersaturation for a 
given crystal-solvent system. 
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An empirical approach also can be used to relate growth 
kinetics to supersaturation by simply fitting growth-rate 
data with a power-law function of the form: 


G = kgo® (25) 


where kg and g are system-dependent constants. Such an 
approach is valid over modest ranges of supersaturation, 
and the power-law function approximates the fundamental 
expressions derived from the above models. 


2. Impurities 


The presence of impurities can alter growth rates substan- 
tially, usually by decreasing them. Furthermore, as de- 
scribed in Section IV.B, impurities can alter crystal mor- 
phology through their effects on the growth rates of crystal 
faces. Mechanisms include: (1) adsorption of an impurity 
on the crystal surface or at specific growth sites such as 
kinks, thereby blocking access to the site by a growth 
unit; (2) formation of complexes between an impurity and 
a growth unit; and (3) incorporation of an impurity into a 
growing crystal and creating defects or repelling the addi- 
tion of a growth unit to the subsequent crystal layer. Few 
of these mechanistic views result in predictive capabili- 
ties, and it is usual to rely on experimental data that are 
often correlated empirically. 

Because impurities most often result in reduced crys- 
tal growth rate, feedstocks to laboratory and bench-scale 
units should be as similar as possible to that expected in the 
full-scale unit. The generation of impurities in upstream 
process units can depend on the way those units are oper- 
ated, and protocols of such units should follow a consistent 
practice. It is equally important to monitor the composi- 
tion of recycle streams so as to detect any accumulation of 
impurities that might lead to a reduction in growth rates. 

The solvent from which a material is crystallized influ- 
ences crystal morphology and growth rate. These effects 
have been attributed to two sets of factors. One has to do 
with the effects of solvent on viscosity, density, and dif- 
fusivity and, therefore, mass transfer. The second factor 
is concerned with the structure of the interface between 
crystal and solvent; a solute—solvent system that has a 
high solubility is likely to produce a rough interface and, 
concomitantly, large crystal growth rates. 


E. Crystal Growth in Mixed Crystallizers 


Population balances on crystals in a crystallizer require a 
definition of growth rates in terms of the rate of change of 
a characteristic dimension: 


dL 


=— (26) 
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Furthermore, the solution of a differential population bal- 
ance requires that the relationship between growth rate 
and size of the growing crystals be known. When all crys- 
tals in the magma grow at a constant and identical rate, the 
crystal—solvent system is said to follow the McCabe AL 
law, while systems that do not are said to exhibit anoma- 
lous growth. 

Two theories have been used to explain growth-rate 
anomalies: size-dependent growth and growth-rate disper- 
sion. As with systems that follow the AL law, anomalous 
growth by crystals in a multicrystal magma produces crys- 
tal populations with characteristic forms. Unfortunately, 
it is difficult to determine the growth mechanism from 
an analysis of these forms. This means that either size- 
dependent growth or growth-rate dispersion may be used 
to correlate population density data without a certainty that 
the correct source of anomalous growth has been identi- 
fied. Determining the actual source of anomalous growth 
is not trivial, but it may be worthwhile since alignment 
between a mathematical model and system behavior en- 
hances the utility of the model. 

Size-dependent crystal growth results when the rate of 
growth depends on the size of the growing crystal. Cer- 
tainly, this may be the case if bulk transport is the control- 
ling resistance to crystal growth, and the literature abounds 
with expressions for the appropriate mass-transfer coef- 
ficients. In the more common situation in which surface 
integration controls growth rate, there are no mechanis- 
tic relationships between growth rate and crystal size, 
and simple empirical expressions are called upon for that 
purpose. 

Growth-rate dispersion is the term used to describe the 
behavior of similar sized crystals in the same population 
exhibiting different growth rates or growth rates that vary 
with time. The consequences of growth-rate dispersion are 
illustrated in Fig. 6, which shows the growth of a crystal 
population that has been immersed in a supersaturated so- 
lution. The spread of the distribution increases as the crys- 
tal population grows; the slower growing crystals form the 
tail of the advancing distribution while the faster growing 
ones form the leading edge. If all crystals in the popula- 
tion grew at the same rate, the distribution would advance 
uniformly along the size axis. Two causes of growth-rate 
dispersion have been observed. In one, the growth rate of 
each crystal in a population is nearly constant, but crystals 
in the population may grow at a different rate; in the other, 
the growth rate of an individual crystal fluctuates about a 
mean value. 

The consequences of anomalous growth depends upon 
the process involved, and this will be pointed out in the 
discussion on population balances. 
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FIGURE 6 Transient population density plot showing growth-rate 
dispersion. 


IV. PURITY, MORPHOLOGY, AND 
SIZE DISTRIBUTIONS 


Crystal properties can be divided into two categories: those 
based on the individual crystal and those involving all crys- 
tals of a given population. The three characteristics of the 
section title compose what is often called crystal qual- 
ity. They represent, along with yield, the most important 
criteria in the operation of a crystallizer. In the following 
discussion, some of the factors that influence purity and 
morphology are described and an introduction is given to 
methods of quantifying crystal size distributions. 


A. Crystal Purity 


The purity of a crystalline product depends on the nature 
of the other species in the mother liquor from which the 
crystals are produced, the physical properties of the mother 
liquor, and the processing that occurs between crystalliza- 
tion and the final product (downstream processing). Im- 
purities can find their way into the final product through 
a number of mechanisms: the formation of occlusions, 
trapping of mother liquor in physical imperfections of the 
crystals or agglomerates, adsorption of species onto crys- 
tal surfaces, as part of chemical complexes (hydrates or 
solvates), or through lattice substitution. 
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Occlusions find their way into the crystal structure when 
the supersaturation in close proximity to the crystal surface 
is high enough to lead to an unstable surface. Such insta- 
bility leads to the creation of dendrites, which then join 
to trap mother liquor in pools of liquid within the crystal. 
Occlusions are often visible and can be avoided through 
careful control of the supersaturation in the crystallizer. 

Mother liquor can be flushed from a cake of crystals 
on a filter or centrifuge by washing with a liquid that also 
may dissolve a small portion of the cake mass. To be ef- 
fective, the wash liquid must be spread uniformly over 
the cake and flow through the porous material without 
significant channeling. Such washing is hindered when 
the crystals themselves have significant cracks, crevices, 
or other manifestations of breakage or the mother liquor 
has a viscosity that is significantly greater than the wash 
liquid. In the latter event, significant channeling (also 
called fingering) may reduce the effectiveness of the wash 
process. 

Lattice substitution requires that the incorporated im- 
purity be of similar size and function to the primary crys- 
tallizing species. In other words, the impurity must fit 
into the lattice without causing significant dislocations. 
An example of such a system is found in the crystalliza- 
tion of L-isoleucine in the presence of trace quantities of 
L-leucine. The two species have similar molecular struc- 
tures, differing only by one carbon atom in the position 
of a methyl side group. In this system, the incorporation 
of L-leucine in L-isoleucine crystals is proportional to the 
concentration of L-leucine in the mother liquor. Moreover, 
the shape of the recovered crystals changes as the content 
of L-leucine in recovered crystal increases. 


B. Crystal Morphology 


Both molecular and macroscopic concepts are important 
in crystal morphology. Molecular structures (i.e., the ar- 
rangements of molecules in specific lattices) can greatly 
influence the properties of a crystalline species and varia- 
tions from a single structure lead to the prospect of poly- 
morphic systems. In such systems, the molecular species 
of the crystal can occupy different locations depending 
on the conditions at which the crystal is formed, and both 
microscopic and macroscopic properties of the crystal can 
vary depending on the polymorph formed. There is, in gen- 
eral, a single stable polymorph for prevailing conditions, 
but that polymorph may not have been formed during the 
crystallization process. In such cases, system thermody- 
namics will tend to force transformation from the unstable 
polymorph to the stable one at rates that may vary from 
being nearly instantaneous to infinitely slow. Additional 
discussion of the molecular structures of crystalline ma- 
terials has been provided elsewhere. 


103 


The characteristic macroscopic shape of a crystal re- 
sults, in large measure, from the internal lattice structure; 
surfaces are parallel to planes formed by the constituent 
units of the crystal. Moreover, although the Law of Con- 
stant Interfacial Angles is a recognition that angles be- 
tween corresponding faces of all crystals of a given sub- 
stance are constant, the faces of individual crystals of that 
substance may exhibit varying degrees of development. 
As a result, the general shape or habit of a crystal may 
vary considerably. 

Crystal morphology (i.e., both form and shape) affects 
crystal appearance; solid—liquid separations such as filtra- 
tion and centrifugation; product-handling characteristics 
such as dust formation, agglomeration, breakage, and 
washing; and product properties such as bulk density, 
dissolution kinetics, catalytic activity, dispersability, and 
caking. 

The shape of a crystal can vary because the relative rates 
of growth of crystal faces can change with system con- 
ditions; faster growing faces become smaller than faces 
that grow more slowly and in the extreme may disappear 
from the crystal altogether. For illustration, consider the 
two-dimensional crystal shown in Fig. 7a and the process 
variables that would cause the habit to be modified to the 
forms shown in Figs. 7b and c. The shape of the crystal 
depends on the ratio of the growth rate of the horizontal 
faces, Gp, to the growth rate of the vertical faces, Gy. For 
the shapes shown in Fig. 7, 


ORON 
GyJ» Gy), Gy). 


Growth rates depend on the presence of impurities, 
system temperature, solvent, mixing, and supersaturation, 
and the importance of each may vary from one crystal face 
to another. Consequently, an alteration in any or all of these 
variables can result in a change of the crystal shape. 

Modeling intermolecular and intramolecular inter- 
actions through molecular mechanics calculations has 
advanced significantly in the past decade, and it has pro- 
vided the basis for prediction of the equilibrium shape of 


{b} fe) 
FIGURE 7 Effect of facial growth rates on crystal shape. 
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a known crystalline species. While not yet uniformly suc- 
cessful on a quantitative basis, the definition and modeling 
of crystal lattice potential energy equations has provided 
an understanding of crystal growth and morphology on the 
molecular level. Derivation of external crystal morphol- 
ogy from internal lattice structures via simulation has been 
proven possible for several organic compounds. Numeri- 
cal minimization techniques, coupled with the appropriate 
valence and nonbonded energy expressions, have enabled 
accurate determination of favorable molecular arrange- 
ments within a wide variety of molecular crystals. 

The shape of crystals obtained as a result of following 
a specific crystallization protocol may be unsatisfactory 
and, as a result, methods for modifying the habit of con- 
siderable interest. The predictive capabilities cited in the 
preceding paragraph are of great utility in such an instance 
as they may be used to determine factors leading to the 
unsatisfactory shape and guide subsequent experiments in 
which a more desirable shape is sought. Inevitably, such a 
search involves extensive laboratory or bench-scale exper- 
iments to determine processing variations that will lead to 
a desired crystal shape. 

As an example of the variations in shape that can be ex- 
hibited by a single crystalline material, consider the forms 
of potassium sulfate shown in Fig. 8. Clearly, the process- 
ing characteristics and particulate properties of the differ- 
ently shaped potassium sulfate crystals will vary. 

The mechanisms and variables affecting crystal shape 
can be categorized as follows: 


FIGURE 8 Shapes of K2SO, crystals. [From Mullin, J. W. (1993). 
“Crystallization,” 3rd ed. Butterworth-Heinemann, London. With 
permission.] 
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1. Intrinsic growth rates 


a. Temperature: The growth rates of individual 
crystal faces depend on temperature, typically 
following an Arrhenius rate law: 


AEg 
G = Go exp (- ) (28) 


RT 


If different crystal faces have different activation 
energies, variation of the temperature at which 
crystallization takes place modifies individual 
growth rates to varying degrees and results in a 
modified crystal shape. 

b. Mixing: The intensity of mixing may determine 
the degree to which bulk mass transfer is involved 
in growth kinetics, and this can influence the 
resulting crystal shape. 

c. Supersaturation: The dependence of growth 
kinetics on supersaturation may vary from one 
crystal face to another. Accordingly, different 
prevailing supersaturations can lead to different 
crystal shapes. 


2. Interfacial behavior 


a. Solvent: Different solvents exhibit different 
interactions with crystal faces and can alter crystal 
shape. A change in solvent also can alter the 
stoichiometry of the crystal (e.g., from a hydrated 
to an anhydrate stoichiometry), which can produce 
crystals with quite different morphology. 

b. Surfactants: Addition of a surfactant to a 
crystallizing system can influence the crystal 
shape in a manner illustrated schematically in 
Fig. 9. Here, surfactant molecules are shown being 
attracted to crystal faces in varying ways; the 
hydrophilic head groups favor the horizontal faces, 
while the hydrophobic tail groups are 
preferentially attracted to the vertical faces. A 
growth unit must displace the surfactant to 
approach a growing crystal face. As hydrophilic 
interactions are typically much stronger than 
hydrophobic ones, the growth unit preferentially 
enters the vertical faces and growth in the 
horizontal direction is favored. 


3. Access to growth site 


a. Blockage by species attracted to growth site: 
Impurities may preferentially locate at a kink or 
other favored growth site and block growth at that 
site. A difference in the character of the kink or 
growth site from one face to another could result 
in modification of the crystal shape. 

b. Species partially fitting into crystal lattice: In 
these instances, an impurity molecule is comprised 
of two parts, one that fits into the crystal lattice 
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Crystal 


FIGURE 9 Attachment of surfactant molecules to crystal surfaces. 


and a second that does not. The part that does not 
fit repulses incoming growth units or causes 
significant interatomic stress because of its 
position in the crystal lattice. If such species are 
purposely added to the crystallizing system to 
modify crystal morphology, they are referred to as 
tailor-made additives. 


C. Crystal Size Distributions 


Most crystallization processes produce particles whose 
sizes cover a range of varying breadth. If the particles 
consist of single crystals, the resulting distribution is a 
crystal size distribution (CSD); on the other hand, if the 
particles consist of agglomerates or other combination of 
multiple crystals, the distribution is a particle size dis- 
tribution. In either case, the distribution is expressed in 
terms of either population (number) or mass. The popu- 


lation distribution relates the number of crystals at each 
size to the size, while the mass distribution expresses how 
mass is distributed over the size range. In the following 
paragraphs, methods for describing and using distribution 
functions will be outlined. 

Size distribution is a major determinant of the proper- 
ties of crystalline products, especially appearance, and to 
downstream processing and handling of crystalline mate- 
rials. Solid-liquid separation by filtration or centrifugation 
can be straightforward with a desired CSD, but it can be 
disastrous when an inappropriate one increases resistance 
to liquor flow through a filter or centrifuge cake. Like- 
wise, CSD affects other downstream processing such as 
the removal of impurities and mother liquor by washing, 
dissolution or reaction of the crystals, and transporting or 
storing crystals. 

Crystal size distributions may be expressed by: (1) his- 
tograms, which are the amount or fraction of mass or 
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TABLE II Sieved KNO3 Crystals from Hypothetical 1-Liter Sample 
Sieve L AM; M(L) F(L) Li AN; N(L) m n 
no., i (um) (g/L) (g/L) (frac) (um) (no./L) (frac) (g/um-L) (no./um.-L) 
1 707 0 32.974 1.000 1611 
2 500 7.296 25.678 0.779 603.5 16 1595 0.0352 0.076 
3 354 11.512 14.166 0.430 427.0 70 1525 0.0789 0.480 
4 240 9.011 5.154 0.156 297.0 163 1362 0.0790 1.430 
5 177 3.145 2.009 0.061 208.5 164 1198 0.0499 2.610 
6 125 1.322 0.687 0.021 151.0 182 1016 0.0254 3.500 
7 88 0.462 0.225 0.007 106.5 181 834 0.0125 4.900 
8 63 0.159 0.066 0.002 T53 175 659 0.0064 7.000 
9 44 0.055 0.011 0.000 53.5 171 488 0.0029 9.000 
10 0 0.011 0.000 0.000 22.0 488 0 0.0002 11.100 
Total 32.974 1611 


number over each increment in size; (2) cumulative distri- 
butions, which are the total or fraction of mass or number 
below (or above) a given size; and (3) density functions, 
which are the derivatives (with respect to size) of cumu- 
lative distributions. These definitions will be illustrated 
by considering a hypothetical potassium nitrate system 
from which a 1l-liter slurry sample has been withdrawn, 
filtered, washed, dried, and sieved to give the results shown 
in Table II. 

The first three and the sixth columns give the sieve data 
and should be read as follows: All of the sieved matter 
passed through the 707-um sieve, and 7.296 g remained 
on the 500-um sieve and had an arithmetic average size of 
603.5 um. Similar descriptions can be given for crystals 
that remained on the other sieves and pan (L = 0). The 
total crystal mass recovered was 39.974 g. A histogram of 
the mass distribution from these data is shown in Fig. 10. 
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Histogram of 
Sieved Crystals 


mass of crystals "° 


on each sieve 


AM; (g) 


The method by which crystals are sized gives either 
number or mass of crystals in a given size range. The sieve 
analysis in the above example gives mass distributions, so 
that the histogram is constructed in terms of crystal mass, 
and a cumulative mass distribution, M(L), can be defined 
as the mass of crystals in the sample passing through the 
sieve of size L. In other words, 


L(i) 
M(L) = 5 AM; (29) 
L=0 
Such calculations give the mass of crystals below size L, 
and the results are shown in column 4 of Table II. Column 
5 gives the cumulative mass fraction distribution: 


M(L) 


F(L) = 
( Miotal 


(30) 


; o : z 
sieve Sizes —* 0 44 63 88 125 177 240 354 500 707 


L(m x 106) 


FIGURE 10 Histogram of size distribution from example. 
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Transforming a mass distribution to a number distri- 
bution, or vice versa, requires a relationship between the 
measured and desired quantities. The mass of a single 
crystal, Merys, is related to crystal size by the volume shape 
factor, kyo) (see Eq. (19)): 


Merys = Pkya L? (31) 


Consequently, the number of crystals on a sieve in the 
example, A N;, can be estimated by dividing the total mass 
on sieve i by the mass of an average crystal on that sieve. If 
the crystals on that sieve are assumed to have a size equal 
to the average of the sieve through which they have passed 
and the one on which they are held, L; = (L,;_; + L;)/2, 
then: 
AM; 


AN; = 73 
PkyoL; 


(32) 
Potassium nitrate crystals have a density of 2.11 x 107! 
g/um?, which allows for the determination of the esti- 
mated crystal numbers on each sieve in Table II. A cumu- 
lative number distribution, N(L), and a cumulative num- 
ber fraction distribution, F(L), can be calculated using 
methods similar to those for calculating M(L) and W(L). 

Mass and population densities are estimated from 
the respective cumulative number and cumulative mass 
distributions: 


AM; 


m(L) = AL, (33) 
z AN; 
n(L) = AL, (34) 
So that if AL; — 0, 
dM(L) e9 
m(L) = => M(L)= | mdL (35) 
0 
and 
nity = Snos f nar (36) 
dL 0 


Equations (33) and (34) are used to obtain the last two 
columns of Table II. 

In the example, all of the results are for the given sam- 
ple size of 1 liter and the quantities estimated have units 
reflecting that basis. This basis volume is arbitrary, but use 
of the calculated quantities requires care in defining this 
basis consistently in corresponding mass and population 
balances. The volume of clear liquor in the sample is an 
alternative, and sometimes more convenient, basis. 

Moments of a distribution provide information that can 
be used to characterize particulate matter. The jth moment 
of the population density function n(L) is defined as: 


mj = [ Lin(L)dL (37) 
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From Eq. (37), it can be demonstrated that the total num- 
ber of crystals, the total length, the total area, and the 
total volume of crystals, all in a unit of sample volume, 
can be evaluated from the zeroth, first, second, and third 
moments of the population density function. Moments of 
the population density function also can be used to esti- 
mate number-weighted, length-weighted, area-weighted, 
and volume- or mass-weighted quantities. These averages 
are calculated from the general expression: 

bage (38) 

Mj 

where j = 0 for a number-weighted average, | for a length- 
weighted average, 2 for an area-weighted average, and 
3 for a volume- or mass-weighted average. 

Crystal size distributions may be characterized usefully 
(though only partially) by a single crystal size and the 
spread of the distribution about that size. For example, 
the dominant crystal size represents the size about which 
the mass in the distribution is clustered. It is defined as 
the size, Lp, at which a unimodal mass density function 
is a maximum, as shown in Fig. 11; in other words, the 
dominant crystal size Lp is found where dm/dL is zero. 
(The data used to construct Fig. 11 are from Table II.) As 
the mass density is related to the population density by: 


m = pk,nL?’ (39) 


the dominant crystal size can be evaluated from the pop- 
ulation density by: 


d(nL?) 
dL o 


The spread of the mass-density function about the dom- 
inant size is the coefficient of variation (c.v.) of the CSD. 


a L=Lp (40) 


fi ; - 
© 100 200 200 400 500 600 TOO 


L (um) 
FIGURE 11 Mass-density function with single mode showing 
dominant size. 
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It is defined by: 
o 
c.v. = Ls (41) 
and estimated from the moments of the distribution: 
amie 1/2 
c.v. = | 7 = 1 (42) 
m4 


This is especially useful for systems that cannot be de- 
scribed by an analytical distribution function. 


V. CRYSTALLIZER CONFIGURATION 
AND OPERATION 


Crystallization equipment can vary in sophistication from 
a simple stirred tank to a complicated multiphase column, 
and the protocol can range in complexity from simply al- 
lowing a vat of liquor to cool to the careful manipulation 
required of batch cyclic operations. In principle, the objec- 
tives of these systems are the same: to produce a product 
meeting specifications on quality at an economical yield. 
This section will examine some of the considerations that 
go into the selection of a crystallizer so as to meet these 
objectives. 

One of the first decisions that must be made is whether 
the crystallizer operation is to be batch or continuous. In 
general, the advantages of each type of operation should 
be weighed in choosing one over the other, but more often 
the decision rests on whether the other parts of the process 
are batch or continuous. If they are batch, then it is likely 
that the crystallizer also should be batch. 

The equipment required for batch crystallization can be 
very simple. For example, some crystalline materials are 
produced by simply allowing a charge of hot liquor to cool. 
After the crystals have formed, the magma is discharged 
through a filter or the liquor may be decanted and the set- 
tled slurry filtered. Very large crystals can be obtained by 
allowing encrustations formed on the walls of these crys- 
tallizers to grow undisturbed; after the system has come 
to equilibrium, the liquor is drained and the crystals are 
removed by scraping them from the surface. 

Batch crystallizers can be used in a campaign to produce 
a particular product and in a second campaign to produce 
another product. Generally, it is not possible to operate 
continuous processes in this way. Batch crystallizers can 
handle viscous or toxic systems more easily than can con- 
tinuous systems, and interruption of batch operations for 
periodic maintenence is less difficult than dealing with in- 
terruptions in continuous processes. The latter factor may 
be especially important in biological processes that require 
frequent sterilization of equipment. Batch crystallizers can 
produce a narrow crystal size distribution, whereas special 
processing features are required to narrow the distribu- 
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tion obtained from a continuous crystallizer. The effects 
of operating variables on crystal size distributions will be 
discussed in Section VI. 

The throughput per unit crystallizer volume is greater 
for a continuous system. Batch units have several operat- 
ing steps in a cycle—charging, heating or cooling, crystal- 
lizing, discharging, and cleaning—and the unit production 
rate is based on the total cycle time, even though the for- 
mation of crystals may occur only during a small portion 
of the cycle. 

It may be easier to operate a continuous system so that 
it reproduces a particular crystal size distribution than it 
is do reproduce crystal characteristics from a batch unit. 
Moreover, the coupling of several transient variables and 
nucleation make it difficult to model and control the op- 
eration of a batch crystallizer. 


A. Relationship of Solubility 
to Mode of Operation 


The driving force for crystal formation can be generated 
through a variety of means, including cooling or heating 
to reduce or increase the system temperature, evaporating 
solvent, evaporative (flash) cooling, inducing a chemical 
reaction (when the reaction product is sparingly soluble, 
the process is called precipitation), adjusting pH, salting 
out through the addition of a nonsolvent, direct-contact 
cooling with a refrigerant, or some other means. All of 
these modes of operation can be implemented in either a 
batch or a continuous process. In addition, two or more of 
the modes may be combined to enhance the product yield. 

As discussed in an earlier section, solubility is intrinsic 
to the solute-solvent system, and the relationship of solu- 
bility to temperature often determines the mode by which a 
crystallizer is operated. Recall for example that the solubil- 
ity of NaCl (see Fig. 2) is essentially independent of tem- 
perature, while NazSOz - 10H20 has a solubility that ex- 
hibits a strong dependence on temperature. Consequently, 
cooling a sodium chloride solution cannot generate sig- 
nificant product yield; solvent evaporation is the primary 
mode of NaCl production. On the other hand, reducing 
the temperature of a saturated solution of sodium sulfate 
generates substantial product and may be used alone or in 
combination with evaporation. 

Cooling crystallizers utilize a heat sink to remove both 
the sensible heat from the feed stream and the heat of crys- 
tallization released or, in some cases absorbed, as crystals 
are formed. The heat sink may be no more than the am- 
bient surroundings of a batch crystallizer, or (as is more 
likely) it may be cooling water or another process stream. 

Evaporative crystallizers generate supersaturation by 
removing solvent from the mixture, thereby increasing 
the solute concentration. They may be operated under vac- 
uum, and in those circumstances it is necessary to have a 
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vacuum pump or ejector as a part of the unit. If the boiling 
point elevation—the increase in boiling temperature due 
to the presence of the solute—is low, mechanical recom- 
pression of the vapor obtained from solvent evaporation 
may be used in some cases to produce a heat source to 
drive the operation. 

Evaporative-cooling crystallizers are fed with a liquor 
whose temperature is such that solvent flashes upon feed 
entry to the crystallizer. They typically are operated under 
vacuum, and flashing of solvent increases the solute con- 
centration in the remaining liquor while simultaneously 
reducing the temperature of the magma. The mode of this 
operation can be reduced to that of a simple cooling crys- 
tallizer by returning condensed solvent to the crystallizer 
body. 

Salting-out crystallization operates through the addition 
of a nonsolvent to the magma in a crystallizer. The selec- 
tion of the nonsolvent is based on the effect of the solvent 
on solubility, cost, properties that affect handling, inter- 
action with product requirements, and ease of recovery. 
Adding a nonsolvent to the system increases the complex- 
ity of the process; it increases the volume required for a 
given residence time and produces a highly nonideal mix- 
ture of solvent, nonsolvent, and solute. 

Melt crystallization operates with heat as a separating 
agent, but a crystalline product is not generated in the 
process. Instead, crystals formed during the operation are 
remelted and the melt is removed as the product. Such op- 
erations are often used to perform the final purification of 
products after prior separation units; for example, the pu- 
rity of an acrylic acid feed may be increased from 99.5 to 
99.9%. Melt crystallizers do not require solids handling 
units nor do they utilize solid—liquid separation equip- 
ment. Finally, in some instances the use of melt crystal- 
lization can eliminate the use of solvents, thereby reducing 
the environmental impact of the process. 


B. Crystallizers 


The basic requirements of a crystallization system are (1) a 
vessel to provide sufficient residence time for crystals to 
grow to a desired size, (2) mixing to provide a uniform 
environment for crystal growth, and (3) a means of gener- 
ating supersaturation. Crystallization equipment is manu- 
factured and sold by several vendors, but some chemical 
companies design their own crystallizers based on exper- 
tise developed within their organizations. Rather than at- 
tempt to describe the variety of special crystallizers that 
can be found in the marketplace, this section will provide 
a brief general survey of types of crystallizers that utilize 
the modes outlined above. 

The forced-circulation crystallizer is a simple unit de- 
signed to provide high heat-transfer coefficients in either 
an evaporative or a cooling mode. Figure 12 shows a 
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FIGURE 12 Schematic diagram of a forced-circulation evapora- 
tive crystallizer. 


schematic diagram of an evaporative forced-circulation 
crystallizer that withdraws a slurry from the crystallizer 
body and pumps it through a heat exchanger. Heat trans- 
ferred to the circulating magma causes evaporation of sol- 
vent as the magma is returned to the crystallizer. This 
type of unit is used to control circulation rates and veloc- 
ities past the heat transfer surfaces, and the configuration 
shown is especially useful in applications requiring high 
rates of evaporation. A calandria that provides heat trans- 
fer through natural convection is an alternative to forced- 
circulation systems. 

Scale formation on the heat exchanger surfaces or at the 
vapor—liquid surface in the crystallizer can cause opera- 
tional problems with evaporative crystallizers. These can 
be overcome by avoiding vaporization or excessive tem- 
peratures within the heat exchanger and by properly in- 
troducing the circulating magma into the crystallizer. For 
example, introducing the circulating magma a sufficient 
distance below the surface of the magma in the crystallizer 
prevents vaporization upon re-entry and forces it to occur 
at a well-mixed zone above the point of re-entry. Alter- 
natively, the magma may be introduced so as to induce a 
swirling motion that dislodges encrustations from the wall 
of the crystallizer at the vapor-liquid interface. 

Figure 13 shows a schematic diagram illustrating the 
configuration of a surface cooling (indirect heat transfer) 
crystallizer. Heat can be transferred to a coolant in an 
external heat exchanger, as shown, or in coils or a jacket 
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FIGURE 13 Schematic diagram of a forced-circulation, surface- 
cooling crystallizer. 
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within the crystallizer. An external cooling surface re- 
quires the use of a circulation pump, but this expense may 
be mitigated by obtaining a higher heat-transfer coeffi- 
cient than would result with the use of coils or a jacketed 
vessel. The rate of heat transfer Q from the circulation 
loop of a cooling crystallizer must be sufficient to reduce 
the temperature of the feed and to remove the heat of 
crystallization of the solute. Assuming that no substantial 
crystallization occurs in the heat exchanger and limiting 
the difference between entering and leaving temperatures 
of the circulating magma (Tin— Tout), SO as to minimize for- 
mation of encrustations, the required magma circulation 
rate Mir, can be determined from the equation: 
; Q 

üa [Cp(Tin E Tout) Jeire cy 
where Cp is the heat capacity of the circulating magma. 
The methods by which Q can be evaluated were discussed 
in Section II. It is not uncommon to limit the decrease in 
magma temperature to about 3 to 5°C; therefore, both the 
circulation rate and heat-transfer surface must be large. 

The feed to cooling crystallizers should be rapidly 
mixed with the magma so as to minimize the occurrence 
of regions of high supersaturation. Such regions lead to 
excessive nucleation, which is detrimental to the crystal 
size distribution. The type of pump used in the circulation 
loop also can lead to degradation of the crystal size dis- 
tribution; an inappropriate pump causes crystal attrition 
through abrasion, fracture, or shear, and most commer- 
cial systems use specially designed axial-flow pumps that 
provide high flow rates and low shear. 

Direct-contact refrigeration can be used if either the 
operating temperature of the crystallizer is low in com- 
parison to the temperature of available cooling water or 
there are severe problems with encrustations. In such an 
operation, a refrigerant is mixed with the crystallizer con- 
tents and vaporized at the magma surface. On vaporizing, 
the refrigerant removes sufficient heat from the magma to 
cool the feed and to remove the heat of crystallization. The 
refrigerant vapor must be compressed, condensed, and re- 
cycled for the process to be economical. Moreover, the 
refrigerant must be insoluble in the liquor to minimize 
losses and product contamination. 

Special devices for classification of crystals may be 
used in some applications. Figure 14 shows a draft-tube- 
baffle (DTB) crystallizer that is designed to provide pref- 
erential removal of both fines and classified product. As 
shown, feed is introduced to the fines circulation line so 
that any nuclei formed upon introduction of the feed can be 
dissolved as the stream flows through the fines-dissolution 
heat exchanger. The contents of the crystallizer are mixed 
by the impeller, which forces the slurry to flow in the in- 
dicated direction. A quiescent zone is formed between the 


Crystallization Processes 


TO VACUUM 

SeoLa * EQUIPMENT 
DATER BAROMETRIC 
CONDENSER 


a 


[to HOTWELL 


STEAM — 
HEATING ELEMENT x : 
— PURGE 

CONDENSATE 


CIRCULATING 
PUMP 


FIGURE 14 Draft-tube-baffle crystallizer. [Courtesy of Signal 
Swenson Division.] 


baffle extending into the chamber and the outside wall of 
the crystallizer. Flow through the quiescent zone can be 
adjusted so that crystals below a certain size (determined 
by settling velocity) are removed in the fines-dissolution 
circuit. In the elutriation leg, crystals below a certain size 
are preferentially swept back into the crystallizer by the 
flow of recycled mother liquor; accordingly, larger crys- 
tals, which have a higher settling velocity, are removed 
preferentially from the system. 

A second major type of crystallizer with special chan- 
neling devices is comprised of those having configurations 
like the Oslo crystallizer shown in Fig. 15. The objective 
of this unit is to form a supersaturated solution by evapora- 
tion in the upper chamber and to have crystal growth in the 
lower (growth) chamber. The use of the downflow pipe in 
the crystallizer provides good mixing in the growth cham- 
ber. As shown, the lower chamber has a varying diameter, 
which can provide some internal classification of crystals. 
The lowest portion of the chamber has the smallest diame- 
ter and can be considered perfectly mixed; as the chamber 
diameter increases, the upward velocity of the slurry de- 
creases and larger crystals tend to settle. In principle, only 
small crystals are supposed to leave the chamber in the 
circulating slurry, to flow through the circulation pump, 
and to enter the upper chamber. As the probability of a 
crystal colliding with the impeller decreases with decreas- 
ing crystal size, the internal classification provided by the 
Oslo crystallizer could provide some control of contact 
nucleation. 

Melt crystallizers can be operated in a variety of ways. 
In one, feed enters the crystallizer and contacts a slurry 
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FIGURE 15 Oslo crystallizer. 


of mother liquor and crystals of the desired product. The 
crystals are forced to move in a specific direction by grav- 
ity or rotating blades. As they flow towards the appropriate 
end of the crystallizer, the crystals encounter a heated re- 
gion and are melted. A portion of the melt is removed 
as product, while the remainder flows countercurrently to 
the crystals, thereby providing some refining and remov- 
ing impure adhering liquid. 

In a second method of operation, the feed material is cir- 
culated through a bank of tubes, each of which has a diam- 
eter of up to about 8 cm. The walls of the tubes are cooled, 
and material crystallizes on them throughout a fixed op- 
erating period. At the end of that period, the remaining 
liquid is sent to a holding tank for further processing, and 
then the tubes are heated slowly to cause partial melting 
of the adhering solids. This step is known as “sweating,” 
and the impure “sweated” liquid produced is removed 
from the crystallizer and held for further processing. Fi- 
nally, the product is obtained by adding additional heat to 
the tubes and melting the remaining adhering solids. The 
actual sequencing of these steps and the reprocessing of 
residual and sweated liquids may be quite complicated. 


VI. POPULATION BALANCES AND 
CRYSTAL SIZE DISTRIBUTIONS 


A balance on the population of crystals in a crystallizer 
can be used to relate process variables to the crystal size 
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distribution of the product or intermediate material. Such 
balances are not independent of those on mass and energy, 
and their solution requires an independent expression for 
nucleation kinetics. 

In formulating a population balance, crystals are as- 
sumed sufficiently numerous for the population distribu- 
tion to be treated as a continuous function. One of the key 
assumptions in the development of a simple population 
balance is that all crystal properties, including mass (or 
volume), surface area, and so forth are defined in terms of 
a single crystal dimension referred to as the characteris- 
tic length. For example, Eq. (19) relates the surface area 
and volume of a single crystal to a characteristic length L. 
In the simple treatment provided here, shape factors are 
taken to be constants. These can be determined by simple 
measurements or estimated if the crystal shape is simple 
and known—for example, for a cube karea = 6 and kyo = 1. 

The beginning point for any balance is the following 
statement: 


input + generation — output — consumption 


= accumulation (44) 


where each of the terms may be expressed as a rate or an 
amount. In a population balance, the number of entities 
(such as crystals) is the balanced quantity and each of 
the terms has dimensions of number of crystals per unit 
time for a differential balance or number of crystals for an 
integral balance. The principles involved in formulating 
a balance are outlined in the following sections, and they 
provide guidance in developing corresponding balances 
for systems whose configurations do not conform to those 
described here. 


A. Perfectly Mixed, Continuous Crystallizers 


The balance equation must be constructed for a control 
volume, which for a perfectly mixed crystallizer may be 
assumed to be the total volume of the crystallizer Vr. 
Then, a balance on the number of crystals in any size 
range (say, Lı to Lz = Lı + AL) must account for crys- 
tals that enter and leave that size range by: (1) convective 
flow, (2) crystal growth, (3) crystal agglomeration, and 
(4) crystal breakage. Agglomeration and breakage can 
be detected through careful inspection of product parti- 
cles, and they can be quite significant in some processes. 
For simplicity, however, they will be assumed negligible 
in the present analysis. The rate of crystal growth G will 
be defined as in Eq. (26); i.e., the rate of change of the 
characteristic crystal dimension L: 


_ aL 
~ dt 
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Then, 
growth rate into the size range = V7(Gn),, (45) 
growth rate out of the size range = V7(Gn)_, (46) 
Lz 
removal rate of crystals in the size range = Vou f ndL 
Lı 
(47) 


Lı 
feed rate of crystals in the size range = V; / nin dL 
Ly 


(48) 
Lı 


f : 0 
accumulation rate in the crystallizer = — nVy dL 
t 


Lı 


(49) 


Substituting the terms from Eqs. (46) through (49) into 
Eq. (44) gives: 


La 
ERA f nin dL 
Lı 


L ð Ly 
= Vr(nG)r, + Vou | ndL+ zf nVrdL (50) 
Lı 


Lı 
Manipulation of this equation leads to 


d(nG) + Vout Vinin 2 ðn (51) 
ðL Vr Vr ot 

Equation (51) may be used as a starting point for the 
analysis of any crystallizer that has a well-mixed active 
volume and for which crystal breakage and agglomera- 
tion can be ignored. As an illustration of how the equation 
can be simplified to fit specific system behavior, suppose 
the feed to the crystallizer is free of crystals and that it is 
operating at steady state. Then, ni, =0 and dn/dt=0. 
Now suppose that the crystal growth is invariant with 
size and time; in other words, assume the system follows 
the McCabe AL law and therefore exhibits neither size- 

dependent growth nor growth-rate dispersion. Then, 


a(nG) an 
aL =i 


Defining a mean residence time tT = Vr/ Vout and applying 
the aforementioned restrictions leads to 


(52) 


dn n 0 53 
Ths (53) 
(t is often referred to as the drawdown time to reflect the 
fact that it is the time required to empty the contents from 
the crystallizer.) Integrating Eq. (53) with the boundary 
condition n =n° at L =Q: 


L 
pee (-=) (54) 


G 
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If the crystallizer has a clear feed, growth is invariant, but 
if the magma volume Vr is allowed to vary, the population 
balance gives: 

ðn a(nG) a(n Vr) Vourn 

u Ce h 
The system model that led to the development of the 
last two equations is referred to as the mixed-suspension, 
mixed-product removal (MSMPR) crystallizer. 

Under steady-state conditions, the rate at which crystals 
are produced by nucleation must be equal to the difference 
in rates at which crystals leave and enter the crystallizer. 
Accordingly, for a clear feed, 


=0 (55) 


[0,6] 1 oe) 
V7B° = Vou f ndL > B° = - [ ndL (56) 
0 T Jo 


For crystallizers following the constraints given above, 
B° = n°G (57) 


For a given set of crystallizer operating conditions, nu- 
cleation and growth rates can be determined by measuring 
the population density of crystals in a sample taken from 
either the well-mixed zone of a crystallizer or the prod- 
uct stream flowing from that zone. Sample analyses are 
correlated with Eqs. (54) and (57), and nucleation and 
growth rates are determined from those correlations. The 
sample must be representative of the crystal population 
in the crystallizer (or leaving the well-mixed unit), and 
experience with such measurements is invaluable in per- 
forming this analysis properly. Figure 16 shows a plot of 
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FIGURE 16 Typical population density plot from perfectly mixed, 
continuous crystallizer. 
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typical population density data obtained from a crystal- 
lizer meeting the stated assumptions. The slope of the plot 
of such data may be used to obtain the growth rate while 
the intercept can be used to estimate nucleation rate. 

Many industrial crystallizers operate in a well-mixed or 
near well-mixed manner, and the equations derived above 
can be used to describe their performance. Also, the sim- 
plicity of the equations describing an MSMPR crystallizer 
make experimental equipment configured so as to meet 
the assumptions leading to Eq. (54) useful in determining 
nucleation and growth kinetics. From a series of runs at 
different operating conditions, correlations of nucleation 
and growth kinetics with appropriate process variables can 
be obtained (see, for example, the discussions of Eqs. (18) 
and (25)). The resulting correlations can then be used to 
guide either crystallizer scale-up or the development of an 
operating strategy for an existing crystallizer. 

It is often very difficult to measure supersaturation, 
especially in systems that have high growth rates. Even 
though the supersaturation in such systems is so small 
that it can be neglected in writing a solute mass balance, 
it is important in setting nucleation and growth rates. In 
such instances it is convenient to substitute growth rate 
for supersaturation by combining Eqs. (18) and (25). This 
gives: 


B° = kyucGi MÅN" (58) 


The constant knuc depends on process variables other than 
supersaturation, magma density, and intensity of mixing; 
these include temperature and presence of impurities. If 
sufficient data are available, these variables may be sepa- 
rated from the constant by adding more terms in a power- 
law correlation. knuc is specific to the operating equipment 
and not transferable from one equipment scale to another. 
The system-specific constants i and j are obtainable from 
experimental data and may be used in scale-up, although 
j may vary considerably with mixing conditions. 

As shown by Eq. (54), growth rate G can be obtained 
from the slope of a plot of the log of population density 
against crystal size; nucleation rate B° can be obtained 
from the same data by using the relationship given by 
Eq. (57), with n° being the intercept of the population 
density plot. Nucleation rates obtained by these proce- 
dures should be checked by comparison with values ob- 
tained from a mass balance (see the later discussion of 
Eq. (66)). 

The perfectly mixed crystallizer of the type described 
in the preceding discussion is highly constrained. Alter- 
ation of the characteristics of crystal size distributions 
produced by such systems can be accomplished only 
by modifications of the nucleation and growth kinetics 
of the system being crystallized. Indeed, examination of 
Eq. (54) shows that once nucleation and growth kinetics 
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are fixed, the crystal size distribution is determined in its 
entirety. In addition, such distributions have the following 


characteristics: 


e Mass density function (from Eq. (39)): 


o73 L 
m = pkan’ L” exp| — — (59) 
Gt 
e Dominant crystal size (from Eq. (40)): 
Lp = 3Gt (60) 
e Moments of n (from Eq. (37)): 
mi = i'n? (Gr) +! (61) 


¢ Total number of crystals per unit volume: 
[o0] 
m= f ndL = mo =n°Gt (62) 
0 
e Total length of crystals per unit volume: 
[o0] 
Lr = I nLdL=m,=n(Gty (63) 
0 
e Total surface area of crystals per unit volume: 


o0 
Ay = kaa | nL? dL = kaeam = 2Kareal” (Gt)? 
0 


(64) 
e Total solids volume per unit volume: 


CO 
Vrs = kval / nL? dL = kyim3 = 6kyan? (Gt) 
0 


(65) 
¢ The coefficient of variation of the mass density 
function (from Eq. (42)) is 50%. 

e The magma density Mr (mass of crystals per unit 
volume of slurry or liquor) is the product of the crystal 
density, the volumetric shape factor, and the third 
moment of the population density function: 


Mr = 6pkyin°(Gt)* (66) 


System conditions often allow for the measurement of 
magma density, and in such cases is should be used as 
a constraint in evaluating nucleation and growth 
kinetics from measured population densities. This 
approach is especially useful in instances of 
uncertainty in the determination of population 
densities from sieving or other particle sizing 
techniques. 


B. Preferential Removal of Crystals 


As indicated above, crystal size distributions produced in 
a perfectly mixed crystallizer are constrained by the na- 
ture of the system. This is because both liquor and solids 
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have the same residence time distributions, and it is the 
crystal residence time distribution that gives the popula- 
tion density function the characteristic exponential form 
in Eq. (54). Nucleation and growth kinetics can influence 
the population density function, but they cannot alter the 
form of the functional dependence of n on L. 

Crystallizers are made more flexible by the introduc- 
tion of selective removal devices that alter the residence 
time distributions of materials flowing from the crys- 
tallizer. Three removal functions—clear-liquor advance, 
classified-fines removal, and classified-product removal— 
and their idealized removal devices will be used here to 
illustrate how design and operating variables can be ma- 
nipulated to alter crystal size distributions. Idealized repre- 
sentations of the three classification devices are illustrated 
in Fig. 17. 

Clear-liquor advance from what is called a double draw- 
off crystallizer is simply the removal of mother liquor 
without simultaneous removal of crystals. The primary ac- 
tion in classified-fines removal is preferential withdrawal 
from the crystallizer of crystals of a size below some spec- 
ified value; this may be coupled with the dissolution of the 
crystals removed as fines and the return of the resulting 
solution to the crystallizer. Classified-product removal is 
carried out to remove preferentially those crystals of a 
size larger than some specified value. In the following 
discussion, the effects of each of these selective removal 
functions on crystal size distributions will be described in 
terms of the population density function n. Only the ideal 
solid—liquid classification devices will be examined. It is 
convenient in the analyses to define flow rates in terms 
of clear liquor. Necessarily, then, the population density 
function is defined on a clear-liquor basis. 


clear liquor 


Vor (R a 1) V 


E solid-liquid 
separator 
perfectly 
mixed 
region 


product 
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Clear-liquor advance is used for two purposes: (1) to 
reduce the quantity of liquor that must be processed by 
the solid—liquid separation equipment (e.g., filter or cen- 
trifuge) that follows the crystallizer, and (2) to separate 
the residence time distributions of crystals and liquor. The 
reduction in liquor flow through the separation equipment 
can allow the use of smaller equipment for a fixed produc- 
tion rate or increased production through fixed equipment. 
Separating the residence time distributions of crystals and 
liquor means that crystals will have an average residence 
time longer than that of the liquor. This should, in princi- 
ple, lead to the production of larger crystals, but because 
the crystallizer is otherwise well mixed, the crystal pop- 
ulation density will have the same form as that for the 
MSMPR crystallizer (Eq. (54)). 

The analysis goes as follows: Let Vin, Vet, and Vout 
represent volumetric flow rates of clear liquor fed to the 
crystallizer, of clear-liquor advance, and of output slurry 
respectively. The population density function is given by 
the expression: 


L 
n=n° exp(- ) (67) 


G Tprod 


where Tproa = Vr/ Vout. Increasing Vcr decreases Vout and 
thereby increases the residence time of the crystals in the 
crystallizer. Unless the increase in magma density results 
in significant increases in nucleation, the utilization of 
clear-liquor advance will produce an increase in the dom- 
inant crystal size. Often the increase is much greater than 
that predicted from theory, and it is suspected that this is 
because the stream being removed as clear liquor actu- 
ally contains varying amounts of fines. If this is the case, 
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FIGURE 17 Schematic representations of idealized removal functions. (a) Clear-liquor advance, (b) classified-fines removal, (c) classified- 


product removal. 
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both clear-liquor advance and fines-removal are enhancing 
crystal size. 

As an idealization of classified-fines removal, assume 
that two streams are withdrawn from the crystallizer, one 
corresponding to the product stream and the other a fines- 
removal stream. Designate the flow rate of the clear so- 
lution in the product stream to be Voy; and set the flow 
rate of the clear solution in the fines-removal stream to be 
(R — 1)Vour. Also, assume that the device used to separate 
fines from the larger crystals functions so that only crystals 
below an arbitrary size Lp are in the fines-removal stream 
and that all crystals below size Lp have an equal probabil- 
ity of being withdrawn as fines. Under these conditions, 
the crystal size distribution is characterized by two mean 
residence times, one for the fines and the other for crystals 
larger than Lr, that are related by the equations: 


= (for L > Lp) (68) 
Vout 
Vr T 
= = — fr L< L 69 
TF RV R (for L < Lr) (69) 


where Vr is the total volume of clear solution in the crys- 
tallizer. 

For systems following invariant growth, the crystal pop- 
ulation density in each size range will decay exponentially 
with the inverse of the product of growth rate and residence 
time. For a continuous distribution, the population densi- 
ties of the classified fines and the product crystals must be 
the same at L = Lp. Accordingly, the population density 
for a crystallizer operating with classified-fines removal 
is given by: 


RL 
n=n° exp- 5 | (for L < Lp) (70) 
Gt 
o (R — 1)Lp L 
n=n expl Gr Jef =| (for L > Lp) 
(71) 


Figure 18 shows how the population density function 
changes with the addition of classified-fines removal. The 
lines drawn are for a hypothetical system, but they illus- 
trate qualitatively what can be demonstrated analytically; 
that is, fines removal increases the dominant crystal size, 
but it also increases the spread of the distribution. 

A simple method for implementation of classified-fines 
removal is to remove slurry from a settling zone in the 
crystallizer. The settling zone can be created by construct- 
ing a baffle that separates the zone from the well-mixed 
portion of the vessel—trecall, for example, the draft-tube- 
baffle crystallizer described in Section V—or, in small- 
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FIGURE 18 Population density plot for product from crystallizer 
with idealized classified-fines removal. 


scale systems, by simply inserting a length of pipe or 
tubing of appropriate diameter into the well-mixed crys- 
tallizer chamber. The separation of crystals in the settling 
zone is based on the dependence of settling velocity on 
crystal size. Crystals entering the settling zone and hav- 
ing a settling velocity greater than the upward velocity of 
the slurry remain in the crystallizer. As the cross-sectional 
area of a settling zone is invariant, the flow rate of slurry 
through the zone determines the cut-size Lp, and it also 
determines the parameter R used in Eqs. (69) through (71). 

In a crystallizer equipped with classified-product re- 
moval, crystals above some coarse size Lc are removed 
at arate Z times the removal rate of smaller crystals. This 
can be accomplished by using an elutriation leg, a hydro- 
cyclone, or a screen to separate larger crystals for removal 
from the system. Using the analysis of classified-fines re- 
moval as a guide, it can be shown that the crystal popula- 
tion density is given by the equations: 


n = n° exp -4 (for L < Lc) (72) 
Gt 
_ ,. Eee ZL) er Le L 
n = n° exp Ge exp a (for L > Lc) 
(73) 


where T is defined as the residence time Vr/ Vout. Figure 19 
shows the effects of classified-product removal on crys- 
tal size distribution; the dominant crystal size is reduced 
and the spread of the distribution becomes narrower. Note 
that it is impossible for crystals smaller than Lc to leave 
the idealized classified-product crystallizer illustrated in 
Fig. 17c. Accordingly, the population densities shown on 
Fig. 19 for the classified-product crystallizer represent 
conditions inside the perfectly mixed region of the unit. 
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FIGURE 19 Population density plot for crystals in crystallizer with 
idealized classified-product removal. 


If both fines and product are removed on a classified ba- 
sis, the population density will be given by the equations: 


n 


RL 
n° exp] - a (for L < Lp) (74) 
T 


z (R — 1)Lf L 
n exo| Gt Jef =| 


(for Lp < L < Lo) (75) 


o exo| (R — a exp| = | 
n=n = — 
Gt Gt 


x 
II 


x epl -2 (for L > Le) (76) 


Selection of a crystallizer that has both classified-fines 
and classified-product removal is done to combine the 
best features of each: increased dominant size and nar- 
rower distribution. Figure 20 illustrates the effects of 
both removal functions on population density. Note that 
this plot of population density results from sampling the 
magma within a crystallizer, not from sampling the prod- 
uct stream, which for the ideal classification devices con- 
sidered here can only have crystals larger than Lc. As 
discussed earlier for the classified-product crystallizer, 
the population densities shown in Fig. 20 represent those 
found in the crystallizer. 

The model of the crystallizer and selective removal de- 
vices that led to Eqs. (74) through (76) is referred to as 
the R-Z crystallizer. It is an obvious idealization of actual 
crystallizers because of the perfect cuts assumed at Lp and 
Lc. However, it is a useful approximation to many systems 
and it allows qualitative analyses of complex operations. 

Although many commercial crystallizers operate with 
some form of selective crystal removal, such devices can 
be difficult to operate because of fouling of heat-exchanger 
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surfaces or blinding of screens. In addition, classified- 
product removal can lead to cycling of the crystal size dis- 
tribution. Often such behavior can be minimized or even 
eliminated by increasing the fines-removal rate. 

Moments of the population density function given by 
Eqs. (74) through (76) can be evaluated in piecewise 
fashion: 


Lr Le oa) 
mi = i LindL +f LindL +f LindL (77) 

0 Lr Le 
Equation (77) is used to estimate the moments of the 
population density function within the crystallizer, not 
of the product distribution. (Recall that moments of the 
distribution within the crystallizer are often required for 
kinetic equations.) Assuming perfect classification, mo- 
ments of the product distribution can be obtained from the 

expression: 


[o0] 
Mi, prod = f LindL (78) 
Lc 
Moments can be used to characterize the material pro- 
duced from or contained in a crystallizer with classified- 
fines or classified-product removal or to evaluate the effect 
of these selective removal functions on product character- 
istics. All that is required is the use of the equations derived 
earlier to relate special properties, such as coefficient of 
variation to the operational parameters R and Z. 


C. Batch Crystallization 


As with continuous crystallizers, the mode by which su- 
persaturation is generated affects the crystal yield and size 
distribution; however, it is the rate at which such super- 
saturation is generated that is most important in determin- 
ing product characteristics. Furthermore, there are infinite 
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FIGURE 20 Population density plot for crystals in crystallizer with 
idealized classifiedfines and classified-product removal. 
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possibilities in selecting cooling profiles, T(t), or vapor 
generation profiles, V(t), or time dependencies of precip- 
itant or nonsolvent addition rates. 

For illustrative purposes, consider that the protocol for 
a cooling crystallizer can involve either natural cooling— 
cooling resulting from exposure of the crystallizer con- 
tents to a heat sink without intervention of a control 
system—or manipulation of cooling to reduce the sys- 
tem temperature in a specific manner. In both cases, the 
instantaneous heat-transfer rate is given by: 


Q = VA(T — Tsink) (79) 


where U is a heat-transfer coefficient, A is the area avail- 
able for heat transfer, T is the temperature of the magma, 
and Tyink is the temperature of the cooling fluid. If Tyink is 
a constant, the maximum heat-transfer rate and, therefore, 
the highest rate at which supersaturation is generated are 
at the beginning of the process. This protocol can lead to 
excessive primary nucleation and the formation of encrus- 
tations on the heat-transfer surfaces. 

The objective of programmed cooling is to control the 
rate at which the magma temperature is reduced so that su- 
persaturation remains constant at some prescribed value, 
usually below the metastable limit associated with pri- 
mary nucleation. Typically the batch is cooled slowly 
at the beginning of the cycle and more rapidly at the 
end. An analysis that supports this approach is presented 
later. In size-optimal cooling, the objective is to vary 
the cooling rate so that the supersaturation in the crys- 
tallizer is adjusted to produce an optimal crystal size 
distribution. 

Protocols similar to those described above for cool- 
ing crystallizers exist for crystallization modes involving 
evaporation of solvent and the rate at which a non solvent 
or a reactant is added to a crystallizer. 

A population balance can be used to follow the devel- 
opment of a crystal size distribution in batch crystallizer, 
but both the mathematics and physical phenomena being 
modeled are more complex than for continuous systems 
at steady state. The balance often utilizes the population 
density defined in terms of the total crystallizer volume, 
rather than on a specific basis: ñ = n Vr. Accordingly, the 
general population balance given by Eq. (51) can be mod- 
ified for a batch crystallizer to give: 


a(n Vr) £ Ə(GnVr) _ dA 2 (GA) | 
or ƏL ðt aL 


The solution to this equation requires both an initial con- 
dition (ñ at t = 0) and a boundary condition (usually ob- 
tained by assuming that crystals are formed at zero size): 


Bt) 
G(0, t) 


0 (80) 


A(0, t) = A° (t) = 


(81) 
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The identification of an initial condition associated with 
the crystal size distribution is very difficult. If the system 
is seeded, the initial condition becomes: 


n(L, 0) = ñseea( L) (82) 


where ñseeq is the population density function of the seed 
crystals. If the system is unseeded, the nuclei often are 
assumed to form at size zero. 

The rate of cooling, or evaporation, or addition of dilu- 
ent required to maintain specified conditions in a batch 
crystallizer often can be determined from a population- 
balance model. Moments of the population density func- 
tion are used in the development of equations relating the 
control variable to time. As defined earlier, the moments 
are 


© ý 
mi =f L'ñdL (83) 
0 
Recognizing that the zeroth moment is the total number 
of crystals in the system, it can be shown that: 


d 
dmi ghia) ge ONT 
dt dt 


Moment transformation of Eq. (80) leads to the following 
relationship: 


(84) 


om j 
— = jGm;_ 85 
at Jmj- ( ) 
Combining Eq. (85) with the relationships of moments 
to distribution properties developed in Section VI.A for 
j=l, 2, 3 gives: 


d mo=Nr dL 

A Gig EGN (86) 
dt dt 

d m=L7 dA 

= = 2Gm, irae 2a = 2G KareaLt (87) 
dt dt 

dm; 


kacama=Ar AM kyo 
= = A 
P7 3Gm, ~— 7 360 ) rt (88) 


Karea 

where Nr is the total number of crystals, Lr is total crystal 
length, Ar is total surface area of the crystals, and My is 
the total mass of crystals in the crystallizer. In addition to a 
population balance, a solute balance must also be satisfied: 


d(VrC) dMy 
dt dt 


where Vy is the total volume of the system, and C is solute 
concentration in the solution. 

The above equations can be applied to any batch crys- 
tallization process, regardless of the mode by which su- 
persaturation is generated. For example, suppose a model 
is needed to guide the operation of a seeded batch crys- 
tallizer so that solvent is evaporated at a rate that gives 


=0 (89) 
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a constant crystal growth rate G and no nucleation; in 
other words, supersaturation is to be held constant and 
only those crystals added at the beginning of the run are in 
the crystallizer. Model development proceeds as follows: 
combining the solute balance, Eq. (89), with Eq. (88), 


d(C), 3p Atk _ 
dt Karea 7 


Recognizing that the process specification requires C to 
be a constant and taking the derivative of Eq. (90): 


0 (90) 


d’ Vr kyo dAr 
C 3 G — =0 91 
g yee (=) di on 
4 Eq. (87) 
eee CoP in= 0 (92) 
dt2 PKyol TL = 
Taking the derivative of the last equation: 
d’ Vr dLr 
C 6pkyiG? —— =0 93 
dt? eee dt 
4 Eq. (86) 
Vr 3 
C PE + 6pky1G° Nr = 0 (94) 


Suppose that the batch crystallizer is seeded with a mass 
of crystals witha uniform size of Leea. The number of seed 
crystals is Nseea, and, as the operation is to be free from 
nucleation, the number of crystals in the system remains 
the same as the number of seed crystals. The initial val- 
ues of total crystal length, total crystal surface area, total 
crystal mass, and system volume are 


Ly(0) = NyceaL seed (95) 
Ar(O) = karea seed L220 (96) 
Mrz(0) = pkvoiNcceaL veg (97) 
Vr(0) = Vro (98) 


On integrating Eq. (94), the following dependence of sys- 
tem volume on time can be obtained: 


C (Vro — Vr) = kvo Nseea[ (Gt)? + 3(Gt)” seca 


seed 


+ 3(Gt)L;.04] (99) 


Therefore, for the specified conditions, the evaporation 
rate (—d Vr/dt) is a parabolic (second-order) function of 
time, and the rate of heat input to the crystallizer must be 
controlled to match the conditions called for by Eq. (99). 

If a cooling mode is used to generate supersaturation, an 
analysis similar to that given above can be used to derive 
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an appropriate dependence of system temperature on time. 
The result depends upon the relationship of solubility to 
temperature. If that relationship is linear, the cooling rate 
varies with time in a parabolic manner; i.e., 


dT F 
-7 = Cut? + Cat + Cs (100) 


An approximation to the temperature-time relationship 
that serves as a good starting point for establishing a fixed 
protocol is given by: 


3 
T=- m- Tad (£) (101) 


where t is the overall batch run time. 

It is clear that stringent control of batch crystallizers 
is critical to obtaining a desired crystal size distribution. 
It is also obvious that the development of a strategy for 
generating supersaturation can be aided by the types of 
modeling illustrated above. However, the initial conditions 
in the models were based on properties of seed crystals 
added to the crystallizer. In operations without seeding, 
initial conditions are determined from a model of primary 
nucleation. 


D. Effects of Anomalous Growth 


Throughout this section, crystals have been assumed to 
grow according to the McCabe AL law. This has simpli- 
fied the analyses of both continuous and batch crystal- 
lizers and, indeed, crystal growth often follows the AL 
law. However, as outlined in Section III, size-dependent 
growth and growth-rate dispersion contribute to devia- 
tions from the models developed here. Both of these phe- 
nomena lead to similar results: In continuous, perfectly 
mixed crystallizers, the simple expression for population 
density given by Eq. (54) is no longer valid. Both size- 
dependent growth and growth-rate dispersion due to the 
existence of a random distribution of growth rates among 
crystals in a magma lead to curvature in plots of Inn vs. 
L. Models for both causes of this behavior exist but are 
considered beyond the scope of the present discussion. In 
batch crystallization, the effects of anomalous growth lead 
to a broadening of the distribution, as was illustrated in 
Fig. 6. 


E. Summary 


The discussion presented here has focused on the princi- 
ples associated with formulating a population balance and 
applying simplifying conditions associated with specific 
crystallizer configurations. The continuous and batch sys- 
tems used as examples were idealized so that the principles 
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could be illustrated, but the concepts can be applied to 
more complicated configurations. Additionally, there has 
been a growing body of work on aspects of population 
balance formulation that greatly extends the ability to de- 
scribe complex systems. Such work has involved anoma- 
lous crystal growth, crystal agglomeration, and crystal 
breakage and necessarily results in substantially more 
complex models. 
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